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ABSTRACT 

Both the detection of galaxies and the derivation of the luminosity function depend upon 
isophotal magnitudes, implicitly in the first case, and explicitly in the latter. However, unlike 
perfect point sources, the fraction of a galaxy's light contained within the limiting isophote is a 
function of redshift, due to the combined effects of the point spread function and cosmological 
dimming. This redshift variation in the measured isophotal luminosity can strongly affect the 
derived luminosity function. Using simulations which include the effects of seeing upon both 
disk and elliptical galaxies, we explore the size of the systematic biases which can result from 
ignoring the redshift variation in the fraction of detected light. We show that the biases lead to 
underestimates in the normalization of the luminosity function, as well as changes in shape. The 
size of the bias depends upon redshift, and thus can mimic galaxy evolution. Surprisingly, these 
biases can be extremely large without affecting (V/V ma x}- However, these biases can be detected 
in the full distribution of V/V max , and in fact may have already been detected in recent surveys. 
Because the systematic biases result from the redshift variation in the fraction of lost light, the 
biases are not significant when the fraction of lost light is always small over the entire survey 
volume, for all galaxy types. However, as modern galaxy surveys now reach higher redshifts, 
lower surface brightnesses, and smaller angular sizes, the effects of seeing and galaxy visibility 
are becoming increasingly important and need to be taken into account. We show that the 
expected biases are not necessarily eliminated when using aperture magnitudes, FOCAS "total" 
magnitudes, or Kron magnitudes, but may be reduced significantly if Petrosian magnitudes are 
used. These considerations may also apply to samples of clusters selected in X-rays. 

1. Introduction 

The luminosity function remains one of the principal tools for quantifying the population of galaxies. 
The evolution of the luminosity function with redshift provides a critical test for assessing changes in the 
galaxy population with time. Variations in the galaxy luminosity function with environment, emission line 
strength, and morphological type all provide broad clues towards distinguishing among galaxy formation 
scenarios. Given that the usefulness of the luminosity function as a cosmological tool depends upon 
comparing luminosity functions at different times or for different galaxy types, it is extremely important 
that the measurement of the luminosity function be unbiased with redshift, galaxy morphology, or survey 
technique. 



1 e-mail address: jd@ociw.edu 
2 Hubble Fellow 



-2 - 



Standard approaches for calculating the luminosity function involve mathematical variants of what 
Binggeli et al. (1988) refer to as the "classical method", wherein each galaxy is weighted by the inverse 
of the maximum volume over which it could have been detected, given the galaxy's luminosity and the 
magnitude limit of the survey. Other methods which reduce the impact of spatial inhomogeneity on the 
derived luminosity function, such as Sandage et al.'s (1979) parametric maximum-likelihood method, 
calculate a likelihood based upon comparing the observed luminosity of each galaxy in the sample to the 
range of potentially observable luminosities at each galaxy's redshift. 

Both of these methods require estimating the range of absolute magnitudes which can be seen at 
each redshift, given the selection criteria of a survey. Any systematic error in this procedure can lead to 
a systematic mismeasurement of the luminosity function. As we will show here, typical procedures for 
measuring the luminosity function do indeed suffer from these systematic biases, largely due to a simplistic 
estimate of how a galaxy's apparent magnitude varies with redshift. 

Most determinations of the luminosity function assume that an isophotal magnitude is a good 
representation of the galaxy's total magnitude, and that the isophotal magnitude varies with redshift as 
the inverse square of the luminosity distance (modulo fc-corrections). While these are perhaps reasonable 
approximations for point sources in deep data with exceptional seeing, or for nearby surveys of extremely 
bright, high surface brightness galaxies, these assumptions break down for real galaxies in many modern 
surveys (e.g. McGaugh 1994). Some of a galaxy's light is lost beyond the outer isophote of a galaxy image, 
due to the combined effects of seeing and surface brightness variations (both cosmological and intrinsic), 
leading to the apparent isophotal magnitude falling off more rapidly with redshift than predicted. Thus, 
field galaxy surveys can easily overestimate the effective volume of the survey, and therefore systematically 
underestimate the luminosity function. More importantly, the degree to which the derived luminosity 
function is underestimated is a strong function of galaxy size and surface brightness, as well as being 
redshift dependent. This can lead to apparent redshift evolution in the galaxy population, and to erroneous 
conclusions about the variation of the morphological mix within the galaxy population. The measured 
luminosity function can also depend upon seeing such that similar surveys with different imaging data 
(e.g. ground based vs. HST) can produce different results. While these biases can be mitigated against 
by using exceptionally deep imaging data or large aperture magnitudes, the majority of deep field surveys 
have used isophotal magnitudes to derive the luminosity function. Instead of isophotal magnitudes, some 
fraction of surveys have used corrected "total" magnitudes which choose an aperture size dynamically, 
based upon either the isophotal area or the intensity- weighted first moment radius. However, in the former 
case, the corrected "total" magnitude actually reduces to using an isophotal magnitude with a fainter 
limiting magnitude, as we will show in ^. Furthermore, even when isophotal magnitudes are not used for 
photometry, the selection of galaxies for a survey depends implicitly on the isophotal magnitude limit of the 
survey images, subjecting the derived luminosity function to many of the biases we explore here. 

The importance of correcting for light lost beyond the limiting isophote has been long known (Humason 
et al. 1956), and has traditionally been studied through curve-of-growth analyses. Sandage (1961) explicitly 
warns against the use of isophotal magnitudes in cosmological tests, pointing out that they are not a 
consistent measure of galaxy flux at all redshifts. Until the recent era of deep redshift surveys, it has 
been sufficient to use growth curves to correct isophotal magnitudes onto the same magnitude system; 
since surveys covered only a limited range in redshift, once galaxy magnitudes were corrected to the same 
isophotal limit, a constant fraction of light was detected for each morphological type over the entire redshift 
range spanned by the data. In contrast, redshift surveys now routinely reach redshifts of 1 (e.g. CFRS, 
AUTOFIB; Lilly et al. 1995a, Ellis et al. 1996), such that even with a single limiting isophote, the fraction 
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of detected light for a galaxy is not the same at all redshifts. Furthermore, deep redshift surveys typically 
have much deeper limiting isophotes for both detection and photometry than do the local redshift surveys 
to which they compare, and typically do not make a correction onto a common isophotal magnitude limit. 
This discrepancy can lead to an apparent increase in the number density of galaxies with increasing redshift, 
as pointed out by McGaugh (1994), and explored in Ferguson & McGaugh (1995) and Phillipps & Driver 
(1995). Yoshii (1993) used similar techniques as will be outlined in this paper to calculate how the use of 
isophotal magnitudes affect the number counts and redshift distributions of faint galaxies. 

In this paper, we investigate how biases resulting from the use of isophotal magnitudes affects the 
luminosity function in particular. We attempt to provide a guide for estimating the degree to which these 
biases may be present in different galaxy surveys. First, in §|^, we outline the method for calculating 
the true volume of a galaxy survey using the fraction of observed light as a function of redshift, which 
we calculate in §|3| as a function of observational parameters and intrinsic galaxy properties. We quantify 
the degree to which this will produce misestimates of the luminosity function, such that one can identify 
the observational regimes where these effects are not important. In §^ we use artificial galaxy catalogs to 
demonstrate how the measured luminosity function varies when different reconstructions of the luminosity 
function are used. In we discuss ways in which the distribution of V/V max can be used to test for the 
presence of luminosity function biases and in we discuss the merits and drawbacks of alternatives to 
isophotal magnitudes. 



Measurements of the field luminosity function begin by calculating an absolute magnitude for each 
galaxy in a redshift survey, using 



where m is the apparent magnitude of an observed galaxy at redshift z, D^(z) is the luminosity distance 
appropriate for the chosen cosmology and k(z) is the fc-correction, which incorporates changes in magnitude 
due to the stretching of the bandpass and the variation in the observed rest-frame wavelength with redshift. 

The relationship between apparent magnitude and absolute magnitude given in equation [l] is only 
correct in the limit where the apparent magnitude m is an accurate measure of the total flux, regardless of 
redshift or morphological type. More typically, however, the apparent magnitude which is used is either an 
isophotal magnitude (e.g. Lilly et al. 1995b, Ellis et al. 1996, Lin et al. 1996), a "total" magnitude measured 
within some multiple of the isophotal area (Small et al. 1997), a modified Kron (1980) magnitude measured 
within an aperture whose size is determined by the first moment radius of the light visible above some 
limiting isophote (Yee et al. 1996, Lin et al. 1997) or, more rarely, an aperture magnitude (Gardner et al. 
1996, Glazebrook et al. 1994); we will restrict ourselves to the case of isophotal magnitudes for the purposes 
of this paper, since isophotal magnitudes are most commonly used, and other magnitude measures reduce 
to this case; we will consider alternative measures of magnitude in 

For an isophotal magnitude measured within a limiting isophote \iu m , some fraction of the light is lost 
outside the outer isophote. The fraction of detected light, f(z), depends both upon intrinsic properties 
of the galaxy (such as the intrinsic central surface brightness /io, the true absolute magnitude M, the 
two-dimensional shape of the galaxy in the absence of seeing (i.e. its light profile)), as well as observational 
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parameters of the survey, such as the limiting isophote /x^ m and the point spread function. 

The fraction of detected light also depends upon the redshift of the observed galaxy. As a galaxy moves 
to higher redshifts, it suffers two effects which rapidly decreases the fraction of light detected within a fixed 
isophotal limit. First, at large enough redshifts the galaxy appears small compared to the point spread 
function, and begins to lose light beyond the limiting isophote, due to the rapid falloff in the point spread 
function with radius. Second, the apparent surface brightness drops off as (1 + z)~ A due to the difference in 
the redshift dependences of the angular diameter and luminosity distances. As the drop in apparent surface 
brightness becomes significant (a factor of 2 at z = 0.2), a larger fraction of a galaxy's light falls below the 
limiting isophote, again increasing the fraction of lost light. The direction of both of these effects is for the 
apparent magnitude to drop off more quickly with distance than predicted by equation |l|. 

If instead, one includes the correction for lost light, equation [j] becomes 



Including the fraction of detected light f{z) has two effects. First, because / < I , the absolute magnitude 
of each galaxy is calculated to be brighter, which compensates for the unobserved light beyond the limiting 
isophote. Second, when equation ^ is inverted to find the maximum redshift at which the galaxy would be 
included in the sample (i.e. z max such that m, so = mum), then the effective luminosity distance is increased, 
and thus a smaller z max is derived than in the absence of the lost light term. This reduces the maximum 
volume V m axi which in turn increases the derived luminosity density (which scales with 1/V max ). 

In luminosity function calculations which reduce the effects of spatial inhomogeneities (Sandage et 
al. 1979, Efstathiou et al. 1988), the derived luminosity function depends upon maximizing the likelihood 
In A = ln^fMi) - In JmZI(z!) H M )d M , summed over all galaxies in the survey. For these derivations, 
including the lost light in equation |^ reduces the range of possible M, which will change the shape of 
the luminosity function. The normalization of the luminosity function in these methods, like the classical 
methods, depends upon the sum of 1/V max> and will also be higher when the effects of lost light are 
included. 



Given equation |2| for the relationship between absolute magnitude and apparent isophotal magnitude 
as a function of redshift, it is in principle straightforward to calculate the luminosity function for a given 
sample. The reconstructions of the luminosity function which are currently in use can readily be applied, 
but with the substitution of the correct values of z maxi or with the correct range of absolute magnitudes at 
a given redshift. 

The difficulty in calculating the luminosity function becomes calculating the form of f(z) for the 
galaxies within a sample. For a galaxy with a fully specified two-dimensional surface brightness profile, one 
can can calculate f{z) at any redshift by first rescaling the central surface brightness by (1 + z)~ A and 
applying the fc-correction, then calculating the angular size using the angular diameter distance, and finally 
convolving the resulting 2-d angular light distribution with the point spread function. Then, f(z) may 
be calculated by integrating the resulting light distribution over the area where the surface brightness is 
greater than the limiting isophote fiu m , or within an area appropriate to any other magnitude measure. 




3. Calculating the Observed Light Fraction f{z) 
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There are several existing papers which calculate the fraction of detected light f(z), particularly in 
the case of large photographic surveys wherein image saturation is important but seeing can be neglected. 
Disney & Phillipps (1983) have calculated f(z) for circular axisymmetric galaxies with both exponential 
and de Vaucouleurs profiles, and used the resulting lost light fraction to calculate the relative V max for 
galaxies of varying surface brightness, at fixed total luminosity^]. Their analysis includes the effects of 
combined isophotal magnitude limits and angular diameter limits. Davies (1990) has expanded upon the 
work of Disney & Phillipps, by generalizing the calculation of the lost light and V max to include combined 
bulge+disk models, and random inclinations. 

The earlier work of Disney & Phillipps has also been extended to higher redshifts by Phillipps et al. 
(1990), who include the effect of cosmological dimming and fc-corrections in calculating both the fraction of 
detected light and V ma x as a function of intrinsic surface brightness and luminosity, for face-on exponential 
disks in the absence of seeing. Using a simple model, they demonstrate how failure to compensate for 
the variation in the effective survey volumes with variation in galaxy surface brightness can lead to 
underestimating the luminosity function for low surface brightness galaxies. 

At moderate redshifts, however, seeing can be as important as surface brightness in modifying the 
fraction of observed light and biasing V ma x- When the half-light radius of a galaxy approaches the size 
of the point spread function (PSF), the observed light profile becomes dominated by the PSF. This can 
becomes significant even at relatively moderate redshifts {z ~ 0.1). In a separate paper (Dalcanton 1998), 
we derive an analytic expression for the observed light fraction f(z) of axisymmetric exponential disks and 
de Vaucouleurs profile ellipticals in the presence of a realistic Moffat (1969) point spread function, which we 
will use for the duration of this paper; Yoshii (1993) has calculated f(z) for the simpler case of a Gaussian 
point spread function. We parameterize the two types of profiles with a characteristic surface brightness 
Ho and a length scale a. For the exponential profile, /irj is the central surface brightness, and a is the 
exponential scale length. For the de Vaucouleurs profile, we chose a and fi such that the de Vaucouleurs 
galaxy has the same half light radius and total flux as an exponential galaxy with the same parameters. 
With this choice, /xq = /j, e — 2.518 and a — r e /1.679, where [i e and r e are the more typical parameterization^ 
of the de Vaucouleurs profile. 

Using the prescription in Dalcanton (1998), and neglecting the /c-correction for generality, we have 
calculated the quantity AM(z) = — 2.51og/(z), which characterizes the degree to which the apparent 
isophotal magnitude m falls off faster than 1/Dl 2 (z). We have plotted AM{z) for different galaxy sizes 
(solid, dashed, and dotted lines) and characteristic surface brightnesses relative to the limiting isophote 
(A/it = fiurn — fJ-o decreasing from left to right), and for variable seeing (increasing from top to bottom), for 
the cases of exponential (Fig. [jja]) and de Vaucouleurs (Fig. |[b]) light profiles. These plots can be used 
to identify the regimes where a survey may suffer from systematic biases in the measured the luminosity 
function; when AM(z) is large over a survey's redshift range, the apparent magnitude of a galaxy can be 
much fainter than would be predicted by eqn. \j\ alone, and thus, a given galaxy drops below some isophotal 
magnitude limit faster than predicted by eqn. fy. 

For both types of galaxy profiles, the curves in Figure [j] show the expected behavior for Atf(z). At low 
redshifts, AM(z) is constant, because the fraction of lost light does not change dramatically with increasing 
distance. When the redshift becomes large enough, however, a given galaxy becomes small with respect to 



3 The relative Vmax as a function of surface brightness, luminosity, and/or scale length is sometimes referred to as the 
"visibility" in the literature (e.g. Disney & Phillipps 1983), since it reflects the relative contribution of different galaxies to a 
given survey. Thus, galaxies which have a smaller Vmax are less likely to appear in a survey, and are thus less "visible" 
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the point-spread function, and, if the redshift approaches 1, the apparent surface brightness drops rapidly 
due to the (1 + z)~ 4 surface brightness dimming. These two effects drop the apparent flux precipitously, 
increasing AM(z). 

At low redshifts, where AM(z) is constant, the apparent magnitude falls off as expected by eqn |l|, 
and thus, as long as the maximum redshift reached by a survey falls within this flat regime, the standard 
reconstruction of the luminosity function is valid (i.e. eqn. [l])^| In contrast, if the limiting magnitude of 
a survey is very faint, then a significant fraction of galaxies will be at large enough redshifts that AM(z) 
is not constant over the survey volume. In this case, the apparent magnitude will be falling faster than 
expected, and thus the luminosity function will be systematically underestimated. In particular and perhaps 
counterintuitively, the luminosity function is more likely to be underestimated in the bright end, since these 
galaxies will be detectable over the largest range in redshift. On the other hand, the redshift at which 
AM(z) becomes large is smallest for intrinsically small galaxies, and for low surface brightnesses galaxies. 
Thus, in deep enough surveys, the luminosity function may be biased in the faint end. The error in the 
derived luminosity function will also depend upon the characteristic redshift of the survey, and thus, the 
luminosity function calculated for the same type of galaxies can be quite different from survey to survey, 
due to variations in the maximum redshift explored by a survey (i.e. mu m ), variations in the seeing, or 
variations in the limiting isophote ixu m . 

Even in the case of optimal seeing, where AM{z) is constant over a large redshift range, there can 
be large effects on the derived luminosity function. As seen if Figure [I], when a galaxy's central surface 
brightness approaches the limiting isophote, there is a global offset in AM(z), which increases as the surface 
brightness decreases. This offset leads low surface brightness galaxies to have systematically underestimated 
luminosities (see discussions by Disney 1976, Disney & Phillipps 1983, for example), even in the absence 
of seeing effects. Thus, only the very nearest low surface brightness galaxies will be bright enough to be 
within the magnitude limit of a survey, and as such, LSBs have much lower detection volumes, even at low 
redshifts where AM(z) is constant; this argument has been framed as the classic "visibility" problem by 
Disney (1976), Disney & Phillipps (1983), Davies (1990), & McGaugh et al. (1995). Furthermore, if surface 
brightness is correlated with absolute magnitude, as seen in clusters of galaxies (e.g. Binggeli et al. 1984), 
then the underestimate of the luminosity becomes increasingly severe at faint absolute magnitudes. This 
correlation would produce a stretching of the luminosity function, which in turn would lead to a smaller 
faint end slope. The effect is nearly independent of seeing, and has been discussed by Phillips et al. (1990). 

Figure [j] also reinforces the value of having very deep, high-resolution imaging data when selecting 
targets for a redshift survey, such that the magnitude limit of the survey is many magnitudes above the 
limiting magnitude of the imaging data. For such cases, the vast majority of galaxies will have limiting 
isophotes for both detection and photometry which are well below their central surface brightness, and as 
Figure |l| shows, such surveys will have far fewer problems with isophotal magnitude biases, as they will 
contain a very large fraction of most galaxies' light out to z ~ 1. 

Finally, to give some appreciation for the values of Afi = fJ,u m — /j,q and a which might be applicable 
for typical surveys, Figure shows /io and a derived from surface photometry for an assortment of 
galaxy samples, chosen to be representative of the range of galaxy morphologies and luminosities (and 
not to be a definitive compilation on galaxy structural parameters). Bright ellipticals tend to have 
21 < fiQ < 18 Bmag/arcsec 2 , while spirals have /iq > 20 B mag/arcsec 2 , down to the current surface 



4 For example, if one were calculating the luminosity function for high surface brightness galaxies selected from HST imaging 
(e.g. upper left in Figure hi), the reconstructed luminosity function would be valid for all redshifts less than 1. 
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brightness limits of major surveys (fio ~ 24 i? mag/arcsec 2 ). Bulges tend to have the highest surface 
brightnesses (15 < /Uo < 21 £? mag/arcsec 2 ), but also tend to be factors of 10 smaller than ellipticals and 
spirals. However, the lack of small spirals with a < 1 kpc is an artifact of only plotting spirals drawn 
from angular diameter limited surveys, which are strongly biased towards identifying galaxies with large 
physical sizes. We also note that, because of fc-corrections and evolutionary effects, the appropriate value 
of can be a strong function of redshift; for example, while elliptical galaxies might have Afj, = 5 in blue 
bandpasses at low redshift, a value of A// = 1 may be more appropriate at redshifts close to 1. There are 
fewer surface photometry measurements available for other band passes, but for reference, the equivalent to 
the Freeman (1970) surface brightness is 20.1 mag/arcsec 2 in r, and ~ 17.6 mag/arcsec 2 in K for galactic 
disks (Courteau 1996, de Jong 1996b). 

To derive A/u, the surface photometry presented in Figure |^ can be compared with the isophotal limits 
for recent field surveys, which we have attempted to summarize in Table [lj. The Table is provided only 
as a rough guide, given that in many cases, the photometric selection criteria are complicated and/or not 
fully specified, and the reader should refer to the original papers for a more complete description. For 
many of the entries we have estimated the limiting surface brightness for detection based upon the given 
photometric limits, or measures of the sky noise, where possible; only the CFRS and Norris surveys have 
adequately tested this limit using realistic artificial galaxy tests. An examination of Table [j] shows that 
the best infrared field surveys to date have only A/i < 1.5 for pure spiral disks, suggesting that the current 
generation of K band luminosity functions reflect the properties of the high surface brightness elliptical 
population alone. 



4. Luminosity Functions 

The above results suggest that the measured luminosity function can depend sensitively upon the 
survey used to select galaxies, and upon the redshifts and morphologies of the galaxies which have been 
selected. In this section we explore in more detail the specific ways in which the luminosity function may 
be affected, by simulating the entire process of measuring a luminosity function, including both sample 
selection and standard reconstructions of the luminosity function. For this exercise, we will neglect spatial 
density inhomogeneities, and use the "classical" l/V max calculation for deriving the luminosity function; 
more sophisticated maximum-likelihood methods should produce identical results in the limit of smooth 
galaxy distributions. 

In general, reconstruction of the luminosity function proceeds as follows: 

Step 1. A sample of galaxies whose apparent magnitude m is above some threshold mum is selected. There 
may also be additional selections for apparent angular size, particularly for nearby photographic surveys. 
Size criteria will be neglected for the following exercise. 

Step 2. Redshifts are measured for the selected galaxies, and, when combined with the apparent magnitudes 
to, are used to derive absolute magnitudes using equation |l|. 

Step 3. The maximum redshift possible for each galaxy is calculated by solving equation [l] for z maxi using 
the derived absolute magnitude from Step 2, and assuming to = mu m . 

Step 4. The volume V within z max is calculated for each galaxy, and then the inverse is summed over all 
galaxies within absolute magnitude bins, to derive the relative number density as a function of absolute 
magnitude. 
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There are two competing effects which can affect the galaxy luminosity function derived with the 
standard methods. First, due to the lost-light beyond the outer isophote, galaxies will have systematically 
fainter derived absolute magnitudes, because Step 2 assumes that the observed magnitude m is a good 
representation of the total flux from a galaxy. Then, when the artificially faint absolute magnitude is used 
to calculate z max in Step 3, one will estimate that the galaxy has a smaller maximum redshift than if all the 
light were included in the measurement of m. This leads one to underpredict the maximum volume V max , 
and thus calculate larger values of the luminosity function in Step 4. One will also calculate larger values of 
V/Vmax for galaxies which are strongly affected by this bias. The second effect works in the opposite sense. 
As shown in §|[ due to the variation in the fraction of detected light f(z) with redshift, a galaxy's apparent 
magnitude can fall off much faster with redshift than predicted by equation [l] alone (i.e. eqn. |). This leads 
one to overpredict the maximum redshift, and thus underpredict the luminosity function. In this case, the 
derived values of V/V max will be smaller than expected. Because these two effects compete with each other, 
it is possible for a sample to have < V/V max >~ 0.5, in spite of the measured luminosity function being 
strongly distorted from the true underlying distribution. 

To demonstrate the size of these effects, we have modelled the case of a moderately deep redshift 
survey in Figure ^, for the specific case of a limiting isophotal magnitude of mu m = 22 within a limiting 
isophote of nu m = 25 mag/arcsec 2 , for galaxies with different characteristic surface brightnesses, and for 
both exponential and de Vaucouleurs light profiles. The analysis assumes that one is calculating the full 
bivariate luminosity function^], where the galaxies are binned not only in luminosity, but also in intrinsic 
surface brightness 

To create a simulated galaxy sample, we have populated a universe with a uniform density of face-on 
galaxies of a fixed characteristic surface brightness, drawn randomly from a Schechter (1976) luminosity 
function <i>(M) with M* = — 21 and ai um = —1.5. Note that the choice of M and /iq fixes the exponential 
scale length, such that variation in M at fixed surface brightness is effectively a variation in size. This 
is a reasonable approach for spiral disks, which span a continuous range of [1q and a, but is somewhat 
artificial for pure elliptical galaxies, which tend to fall along a locus in fiQ vs a (Figure [2]). We retain this 
choice for ellipticals, however, because it more clearly elucidates the general trends which are expected, 
and because, in general, a bivariate luminosity function will probably not further subdivide a luminosity 
function by profile shape, in addition to surface brightness. For each surface brightness, we calculate 
the observed isophotal magnitude using Moffat convolved galaxy profiles and generate a catalog of 5,000 
galaxies which fall above the isophotal magnitude limit. The resulting catalog is then analyzed using the 
methods most typically applied to galaxy catalogs (Steps 1-4 above; eqn. [I]), and then reanalyzed using the 
proper correction for lost light (eqn. ^|) . For the sake of generality, we have neglected the fc-correction term 
in eqn. ^ both in generating the catalog and in analyzing it; we will discuss the role of the fc-correction 
below. The final luminosity function and V/V max distributions are shown in Figure]^. 

There are several obvious trends demonstrated in the luminosity functions plotted in Figure |3|: 

1. As expected from Figure [l], the errors in the derived luminosity function are greater for low surface 
brightness galaxies, due to the increased fraction of lost light, and the lower redshift at which AM(z) 
becomes large. The biases are slight when the characteristic surface brightness fiQ is 5 mag/arcsec 2 above 



5 This has also been referred to at the bivariate brightness function (BBF) by Phillipps et al. (1990) and Boyce & Phillipps 
(1995). We prefer to retain the phrase "luminosity function", to maintain the analogy to previous work on the galaxy luminosity 
function. See Boyce & Phillipps (1995) for a good discussion of the elements necessary to a a survey designed to measure the 
full bivariate luminosity function 
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the limiting magnitude, significant at A/j, — 3 mag/arcsec 2 , and severe at A/i = lmag/arcsec 2 (where 
A^i = m im - no). 

2. The underestimate in the luminosity function is worse at high redshift, due to the rapid variation in 
the lost light fraction at increasing redshift. Furthermore, in low- redshift subsamples, some portion of the 
sample has V max fixed by the upper cutoff in redshift, which reduces the impact of ignoring the lost light 
corrections. This suggests that when the luminosity function is derived in intervals of redshift, evolution 
detected in the luminosity function may be correct for the lower redshift intervals, but evolution between 
the penultimate and the highest redshift interval is suspect. 

3. The underestimate of the luminosity function is worse at bright magnitudes, as these galaxies tend 
to be at larger redshifts, where the light lost fraction changes rapidly. 

4. The luminosity functions appear to shift to fainter Af* with decreasing surface brightness, due to 
the increasing offset in AM(z) (see Figure [l]). 

5. The mismeasurement of the luminosity function can mimic luminosity and density evolution. The 
luminosity function appears to shift towards brighter magnitudes and to larger normalizations as the 
redshift decreases (see the comparison between the dotted and solid curves in the lower left panels of Figure 
^). Including fc-corrections and evolutionary corrections can induce the opposite effect. Most galaxies are 
probably intrinsically brighter at higher redshifts, and thus the appropriate value of A/^o may be larger 
with increasing redshift, if the evolutionary corrections are larger than the fc-corrections (as might be the 
case for galaxies with an exponentially declining star-formation rate). The luminosity function will then 
shift to higher luminosities at larger redshifts due solely to biases, and will therefore have the appearance 
of more rapid luminosity evolution than is the case. 

6. For samples with the same number of galaxies, the underlying luminosity function has a much 
higher normalization for low surface brightness galaxies, due to their smaller mean distance in isophotal 
magnitude limited surveys; note the difference in the mean redshift as the surface brightness is decreased. 
This is the classical "visibility" effect discussed by Disney (1976), Disney & Phillips (1983), Davies (1980), 
McGaugh et al. (1995), and others. The relative normalizations of the luminosity functions suggest that 
galaxies with A[i = 1 mag/arcsec 2 will be underrepresented by a factor of 10, compared to a population of 
galaxies with an identical space density and Afi = 5 mag/arcsec 2 . If the lower surface brightness galaxies 
are not treated separately, then the derived luminosity function will be 10% higher than would be correct 
for A/x = 5 mag/arcsec 2 , but will be a factor of 10 smaller than is correct for the entire surface brightness 
range spanned by the sample (A/x > lmag/arcsec 2 ). In such instances, one should quote that the derived 
luminosity function is representative of only the high surface brightness fraction of the sample, even if the 
survey is capable of reliably detecting galaxies at much lower surface brightnesses. 

7. At the same surface brightness, the mismeasurement of the luminosity function is comparable for 
both exponential and de Vaucouleurs profiles, as would be expected by the general similarity of the detected 
light fraction f(z) plotted in Figure [D. However, ellipticals and bulges have somewhat higher surface 
brightnesses than spiral disks (Figure ||), by about 1-2 mag/arcsec 2 at the bright end, and thus ellipticals 
and early type spirals will be tend to be less affected by luminosity function biases, when the fc-correction 
is negligible. In contrast, Figure || also shows that ellipticals and spirals have comparable half-light radii, 
and thus both types of galaxies will begin to be affected by seeing at similar redshifts. Bulges have much 
smaller half-light radii than their spiral hosts, and while a high surface brightness bulge will help a spiral 
galaxy's detectability, at redshifts z > 0.1 the bulge profiles begin to be strongly affected by seeing for 
typical ground based observations. As a result, bulges may tend to become undetectable at comparable 
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redshifts as their associated disks (e.g. compare the tracks in Figure [j|b] for Afi = 7 and a = 0.25 kpc to 
the tracks in Figure |[a] for A/i = 3 and a = 3 kpc). 

8. There is less appearance of "evolution" for the de Vaucouleurs profiles than for the exponentials, in 
the absence of the fc-correction, because f(z) is constant over a slightly larger range in redshift for the de 
Vaucouleurs profile, while having a sharper increase at high redshift, and because f(z) varies more strongly 
with surface brightness in exponential disks than in de Vaucouleurs ellipticals. These effects also lead to 
ellipticals having smaller mean redshifts than comparable spirals at high surface brightnesses (A/i > 3) and 
larger mean redshifts at low surface brightnesses (A/i < 1). 



5. V/V max Distributions 

The degree to which the luminosity functions seen in Figure [| have systematic problems can be 
partially diagnosed using the derived V/V max distributions. As presented in Figure^, as the corrections for 
lost light become larger (i.e. at lower central surface brightnesses and higher redshifts) the distributions 
of V/V m ax become sloped, such that the distributions are elevated at small values of V/V max and are 
depressed at larger values. The nature of the deformation is such that it nearly preserves the expectation 
value of < V/V m ax >= 0.5, even in the presence of severe problems in the derivation of the luminosity 
function. Thus, the full distribution of V/V max is needed to test for the presence of lost-light biases. In 
the presence of the necessary /c-correction, however, the full V/V max distribution may well deviate from the 
relatively straight line seen in Figure |3|, and pull the mean < V/V max > from its expected value of 0.5. 

Ideally, the V/Vmax test should be performed in several redshift bins, given that the low redshift 
subsample of the A/i = 1 mag/arcsec 2 sample is clearly biased in its luminosity function derivation, while 
having a nearly normal V/V max distribution. This apparent normality results from V ma x being fixed by the 
upper cutoff in redshift for a large fraction of the galaxies. In general, the slope in the V/V max distribution 
is somewhat more severe for disk galaxies than comparable elliptical galaxies. 

We note that the expected signature has already been seen in the data used for the Autofib survey 
(Ellis et al. 1996; their Figure 2). Even after correction for redshift incompleteness, their published V/V m ax 
distributions systematically fall off towards large values of V/V maxi and look remarkably similar to the 
plots in Figure 0. Most other existing surveys have published only the mean < V/V max > and cannot be 
immediately tested for the biases discussed here. 



6. Alternative Measures of Luminosity 

The systematic biases in the luminosity functions presented in Figure |^ are a direct result of the 
practice of using isophotal magnitudes to measure galaxies' luminosities. There are other ways to measure 
galaxy magnitudes however, which, while less commonly used, can potentially improve the accuracy of the 
luminosity function reconstruction. 

The easiest alternative to the isophotal magnitude is the aperture magnitude. One can choose a large 
enough aperture to safely contain a large fraction of a typical galaxy's light, and thus have a measure of the 
total luminosity which is not biased with surface brightness. However, aperture magnitudes are strongly 
biased with galaxy size; physically large galaxies will have a smaller fraction of their light detected than 
will smaller galaxies, leading the brighter galaxies to be measured with a smaller intrinsic luminosity. The 
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degree of mismeasurement will also be a function of redshift, with a larger fraction of a galaxies' light being 
detected at higher redshift. This redshift dependence is opposite to the behavior of isophotal magnitudes, as 
can be seen in the top row of Figure ^, where we have plotted the error in the measured absolute magnitude 
as a function of redshift, for a 3" diameter aperture (as in Lilly et al. 1991, for 1" seeing), applied to various 
galaxy profiles and physical sizes. Because it is the redshift dependence of AM which led to the artificial 
appearance of galaxy evolution in Figure ^, aperture magnitudes should also produce artificial evolution, 
although with the opposite luminosity evolution. As can be seen from Figure |5|, even for a large aperture 
(3xFWHM), a typical spiral with a = 3kpc is biased by over 1 magnitude a z = 0.1. 

Another common practice is to attempt to correct a measured isophotal magnitude to a "total" 
magnitude. The "total" magnitudes provided by the FOCAS package (Valdes 1982) are created by 
extending the isophotal area by a factor of 2, and then measuring the flux within the larger aperture. 
While not immediately apparent, the FOCAS "total" magnitudes are actually equivalent to an isophotal 
magnitude measured at a fainter limiting isophote. Thus, one can consider using a FOCAS "total" 
magnitude as increasing the value of A/i = fiu m — /io, where \iu m is the limiting isophote for the 
initial detection of the galaxy. For pure exponentials, the new effective magnitude difference is always 
A/x e // = VAA/i, where A is the factor by which the isophotal area is grown. For pure de Vaucouleurs 
profiles, A/x e // = A 1 / 8 A/i + 5.81 (A 1 / 8 — 1), and for galaxies which are small enough to have nearly Gaussian 
profiles due to the seeing, A/z e // = AA/.1 — 5(A — 1) log er/a, where a is the standard deviation of the 
Gaussian point spread function, and a is the galaxy's scale length. Because the increase in A/j, e ff is largely 
a multiple of the isophotal value of A/i, changing from an isophotal magnitude to a "total" magnitude does 
not greatly reduce the level of systematic bias for lower surface brightness galaxies; a galaxy with A/i =1.5 
has A/j, e ff = 1.5 — 2.1 when FOCAS "total" magnitudes are used. 

There are other forms of "total" magnitudes, besides those used by FOCAS. One of the most common 
of these is the Kron (1980) magnitude, which sets the aperture to be a multiple of the intensity- weighted 
first moment radius r\ , instead of a multiple of the isophotal area. With careful consideration of the effects 
of seeing, the Kron magnitude can provide a fairly robust method of providing a magnitude measurement 
which is consistent over a range in redshifts and observing conditions (see Bershady et al. 1994 for a 
particularly cogent discussion of the relevant issues) . The consistency of the Kron magnitude with variations 
in surface brightness has not yet been sufficiently explored, however. Galaxies whose apparent central 
surface brightness approaches the isophotal limit are typically detected out to only one scale length, and 
thus will tend to give systematically smaller values of r\ ; the degree of this mismeasurement will depend 
sensitively on the details of how r\ is measured (iteratively, directly, etc) . 

We test the sensitivity of Kron magnitudes to variations in surface brightness in Figure ^, where we 
have assumed that T\ is measured above some limiting isophote /i/j m , and that the Kron magnitude is 
calculated within a radius r ap — 2ri, as suggested in Kron (1980). Figure f| shows that the Kron magnitude 
underestimates the total luminosity at low surface brightnesses (^u m — Mo 3), much like isophotal 
magnitudes do. For A[i ~ 3, errors in the absolute magnitude are larger than 0.3 — 0.7 m , and increase 
in size towards fainter surface brightnesses (i.e. smaller A/i). The errors at low surface brightnesses are 
also larger when the PSF dominates the galaxy profile. The curves shown in Figure |^ are actually worst 
case scenarios, as we do not calculate r\ iteratively, nor do we prevent the Kron radius 2r x from becoming 
smaller than either some multiple of the original limiting isophotal radius or some fixed aperture radius, 
both of which will give better estimates of the total magnitude than given in Figure f|. However, as in the 
case of FOCAS total magnitudes, calculating r\ iteratively is equivalent to using a fainter value of fium, 
and thus a somewhat larger value of A/i. The curves in Figure will therefore be identical when iterative 
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measures of r\ are used, but the appropriate value of A/i will be larger. The situation is similar when there 
is a minimum allowed Kron radius, if the minimum radius is set to be a multiple of the original limiting 
isophotal radius (as in SExtractor, Bertin & Arnouts 1996). If the minimum Kron radius is a fixed aperture, 
then the biases for the faintest galaxies will become those suffered by aperture magnitudes. 

Another alternative to either the isophotal or aperture magnitude is the Petrosian magnitude (1976). 
To define a Petrosian magnitude, one fixes the value of the parameter 77, finds the Petrosian radius where 
77 =< S(r) > /S(r), and then measures the flux within that radius (here, < S(r) > is the mean surface 
brightness within r). Because the characteristic surface brightness cancels in the definition of 77, the 
Petrosian radius will be identical for galaxies with the same profile shape, independent of their surface 
brightness, and thus the fraction of detected light f(r) will be constant with surface brightness as well. 
Furthermore, unlike for aperture magnitudes, the fraction of a galaxy's luminosity contained within the 
Petrosian radius is independent of a galaxy's size or distance, and depends only upon the particular choice 
of 77, and the shape of the galaxy's profile (in the absence of seeing). These qualities make the Petrosian 
magnitude a nearly ideal choice for measuring luminosity functions in a way that is unbiased with redshift 
or surface brightness. The stability of the Petrosian magnitude can be seen in the bottom of Figure ||. 
The Petrosian absolute magnitude is nearly constant with redshift for pure exponential disks of different 
sizes, even when a 1" PSF is included. For de Vaucouleurs profiles, the absolute magnitude varies by only 
0.2 magnitudes out to a redshift of 2, and differs from exponential profiles by roughly the same amount. 
We have chosen the value of 77 carefully, to give a good compromise between including a large fraction of 
a galaxy's light, yielding consistent results for elliptical and spiral profiles, and producing small enough 
Petrosian radii that S(r) can be accurately measured; for different choices of 77, the results in Figure [| could 
be very different. The comparison of exponential profiles to de Vaucouleurs profiles in Figure]^ suggests that 
even for a realistic mixture of galaxies, Petrosian magnitudes should be able to make consistent measures of 
luminosity over a large range in galaxy size, surface brightness, and redshift, leading to a vast improvement 
over the standard use of isophotal magnitudes. 

While Petrosian magnitudes appear to offer a nearly ideal tool for measuring galaxy luminosities, 
the plots in Figure || disguise several difficulties. In general, Petrosian magnitudes are subject to large 
errors at low signal-to-noise ratios and for poorly sampled profiles. To find the radius corresponding to a 
particular value of 77 in a faint galaxy image, one must carefully interpolate over a relatively small number 
of pixels to find r, S(r) and < S(r) >. For low surface brightness galaxies or normal galaxies at moderate 
redshifts, the particular value of 77 may require accurately measuring extremely faint surface brightnesses 
(i.e. no ~ 30 mag/arcsec 2 ), which may be well beyond the capabilities of the available imaging data, due 
to sky noise, flat-fielding, crowding, etc. Thus, while Figure [| suggests that Petrosian magnitudes are 
unbiased with surface brightness or size in ideal data, for noisy, binned data, there should be both noise and 
systematic biases which will get worse with decreasing surface brightness and redshift, although hopefully 
to a lesser degree than isophotal magnitudes. 

Wirth (1996) has explored the systematic errors in Petrosian magnitudes for different interpolation 
schemes, and for different signal-to-noise ratios. He finds that the simplest interpolation methods are highly 
biased with signal-to-noise, and while more sophisticated ones can be constructed, they are only useful for 
galaxies which have a signal-to-noise ratio above 10 within their true half-light radii. At SNR0.5 = 10, the 
best estimators have a typical scatter in the Petrosian radius of more than 50%, and a systematic bias of 
about 30%, which will translate into systematic errors and increased scatter in the measured Petrosian 
magnitudes for faint galaxies. The simplest interpolation schemes produce much larger systematic biases, 
which reach over 80% in the Petrosian radius at SNR0.5 = 10. Because of these systematics, the robustness 
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of any adopted scheme for measuring Petrosian magnitudes should be thoroughly tested as a function of 
surface brightness, size, and signal-to-noise, before they can be safely invoked as a panacea for the problems 
which can plague isophotal magnitude determinations of the luminosity function. 

Wirth (1996) also finds that for galaxies with high signal-to-noise (SNR0.5 ~ 80), the measured 
Petrosian radii are good to within 10% which suggests that one could perform a nearly ideal luminosity 
function reconstruction for deep enough imaging data. However, even the best data is subject to the strong 
isophotal magnitude biases discussed in the rest of this paper. The detection of galaxies is fundamentally 
an isophotal signal-to-noise limited process; one will never measure photometry for a galaxy which does not 
have substantial light above the sky noise of a given image. Thus, all surveys have an implicit isophotal 
magnitude limit, which is roughly the standard completeness limit of a survey. In determining the biases in 
a luminosity function measurement, even if using Petrosian magnitudes, one must also consider the imact 
of the implicit isophotal magnitude limit as well. 

7. Conclusions and Complications 

In the current epoch of deep CCD imaging, some surveys have extended into a redshift regime where 
the light lost beyond the outer detection isophote, due to intrinsic surface brightness and seeing, can have 
large, systematic effects upon the derived luminosity function, without violating the standard < V/V max > 
tests which are often used to diagnose incompleteness and unusual systematics. We have attempted to 
give a general guide to determining the regimes where corrections for lost light become important, and to 
estimating the degree to which neglecting the corrections can affect the resulting luminosity function. We 
have also shown how these effects can easily mimic evolution in the galaxy population. As an aside, we 
note that the biases we discuss are potentially relevant to any work which involves detection of extended 
objects, such as the detection of galaxy clusters in X-ray data. 

In our investigation, we have taken the simplest, most general assumptions, in order to clarify the 
general biases which affect luminosity function measurements, and to estimate the size of these biases in 
different regimes. This approach comes at the expense of glossing over some subtleties which can affect real 
surveys. Perhaps the most glaring complication is that in many surveys, the limiting isophote for detection 
is substantially brighter than the limiting isophote for photometry. For example, the CFRS survey (Lilly et 
al. 1995a) uses a very deep isophotal limit of roughly [ii AB = 28.5 for photometry (which guarantees that 
they measure a consistently large fraction of a detected galaxy's light), but has an effective isophotal limit 
for detection of roughly fii AB ~ 25.5 (estimated from the central surface brightness detection limit in their 
Figure 2), due to sky noise and flat-fielding uncertainties. In such cases, one must take additional care in 
calculating V max , such that the maximum redshift is derived with the additional criteria that the object 
can be detected above the limiting isophote for detection, as well as having an isophotal magnitude brighter 
than the limiting magnitude of the survey. Testing for a galaxy's detectability is further complicated by 
the usual practice of smoothing an image before searching for galaxies within it. For imaging data which is 
substantially deeper than the limiting magnitude of the spectroscopic subsample, many of these problems 
can be avoided. 

For simplicity, we have also neglected the mix of complicated morphologies and viewing angles present 
in the real universe, and instead used face-on circularly symmetric exponential disks and dc Vaucoulcurs 
profile ellipticals for illustration. As almost all galaxies can be reasonably approximated by a sum of the 
two profile types, reality is likely to be bracketed by these two limiting cases. If a survey is in a redshift 
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regime or a surface brightness regime where the corrections for lost light are probably large (see Figure 
|l|), then the optimal method for deriving V max is to first fit an inclined bulge/disk model to each galaxy, 
then to calculate V max by simulating the individual galaxies as they would be observed at higher redshifts 
(including the likely fc-corrections), using either analytic calculations or artificial galaxy tests. In the latter 
case, one can readily incorporate the complications of galaxy detection as a function of redshift, as well as 
the effects upon the photometry. We have also shown how the use of Petrosian magnitudes can greatly 
reduce the size of the needed corrections, but have argued that their benefits are most substantial at 
relatively high signal-to-noise ratios. 

While making these further corrections can be a fair bit of work, and will certainly be subject to 
their own errors and biases, for the typical size of galaxy surveys (< 1000 galaxies) such an approach is 
not necessarily untenable, and for larger galaxy surveys, such as the upcoming 2DF and Sloan surveys, 
this can be approached in a binned, statistical way, and/or can be heavily automated. Many groups have 
already implemented routines to estimate structural parameters for galaxies using deconvolution and/or 
maximum likelihood methods (e.g. Schade et al. 1996, Ratnatunga et al. 1994), and from there, it is a small 
extrapolation to calculate the correct V max . At the least, it is preferable to make a first-order correction 
when needed, than to ignore effects which may be large enough to obscure the cosmology that one is hoping 
to illuminate. 
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Table 1. Photometric Parameters for Some Recent Faint Field Surveys 



Survey 


Band 


Seeing 

(") 


Magnitude Limits 


Limiting Surface Brightness 
Photometry Detection 


CFRS a 


Iab 


0.6" - 1" 


17.5 < Iab < 22.5 


28.0 Iab/O" 


25.5 Iab/O" 


LCRS b 


Gunn-r 


~ 1.3" 


15 < r iso < 17.7 


~ 23r/n" 


~ 22r/D" 


CNOC c 


Gnnn-r, Bab 


? 


r < 20.5, 21.5, 22.0 


j 


? 


AutoFib d ' c 


bj 


? 


19.5 < bj < 22 


26.5 bj/n" 


~ 25 6.//D" 


Norris f 




~ 2.0" 


r core < 21.7, r tot < 20 


k 


24.1 r/n" 


Gardner 8 


K l 


? 


K < 15 


1 


22.5 J, 18. 5 K/n" 


Glazcbrook h 


K 


1.9" -3.5" 


K < 17 


in 


~ I9if/n" 



"Lilly et al. 1995a 
b Shectman et al. 1996 
°Yec et al. 1996, Lin et al. 1997 
d Ellis et al. 1996 

°Paramctcrs refer only to AutoFib-faint sample; luminosity function includes data from other surveys with 
different selection criteria 

f Small et al. 1997 

g Gardner et al. 1996 

h Glazebrook et al. 1994 

1 K photometry performed only for galaxies detected in / 
J Variable aperture magnitudes used 
k FOCAS total magnitudes used 
l 10" apertures used 

m 4" — 8" apertures used, depending on seeing. 
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9. Figure Captions 

Fig. 1. — [a] The fraction of the light detected above a limiting isophote, represented as the error in the 
derived absolute magnitude (i.e. AM(z) = — 2.5 log/(z)), is plotted as a function of redshift, for different 
seeing (from top to bottom: FWHM=0.1", 1.0", 2.0", for a Moffat profile with (3 = 5), different galaxy 
surface brightnesses relative to the limiting isophote (left to right: A/i = fiu m — (1q — 7, 5, 3, 1), and different 
galaxy sizes (a — 0.25 kpc, 1 kpc, 3 kpc: dashed, dotted, and solid lines, respectively), for circularly symmetric 
galaxies with exponential surface brightness distributions. f(z) includes the effect of cosmological dimming 
on reducing the fraction of light seen at high redshift, but does not include the /c-correction, which can modify 
the effective central surface brightness at high redshift. We have assumed Ho — 50km/s/Mpc and flo = 1. 
For larger values of the Hubble constant, the upturn in AM(z) will occur at larger redshifts; however, the 
characteristic size of the galaxies in Figure || will become smaller, [b] Same as [a], but for a de Vaucouleurs 
profile. See Figure ^| for a treatment of aperture magnitudes and Petrosian magnitudes. 

Fig. 2. — The values of no and a for a variety of galaxy samples. The solid circles and open squares are 
spiral bulges and spiral disks, respectively, drawn from de Jong (1996a), and the asterisks are ellipticals in 
A157 and A3574 from Jorgensen et al. (1995). The solid triangles are low surface brightness galaxies (LSBs) 
from the UGC catalog from Knezek (1993). The crosses are spiral disks from an assortment of other sources 
(Romanishin et al. 1973, Boroson 1981, van der Kruit 1987, Sprayberry et al. 1995, de Blok et al. 1995, 
McGaugh & Bothun 1994) and are primarily selected from the NGC, UGC, or POSS-II catalog of LSBs 
(Schombert et al. 1992), all of which arc cither explicitly or implicitly angular diameter limited field surveys. 
The stars are the local group dwarf spheroidals compiled by Caldwell et al. 1992. The dashed line is the 
locus of Mb = —20. For the de Vaucouleurs profile galaxies (asterisks and filled circles), the conversion 
from effective radius r e and effective surface brightness fi e to fj,Q and a is as given in the text. All data has 
been converted to the Johnson B magnitude system. The data are chosen to represent the existing range 
of galaxy morphologies and luminosities, and not to argue for a particular relation between a and /Iq for a 
given profile type. 

Fig. 3. — The luminosity functions (left column) and V/V max distributions (right column) for samples of 
galaxies with different central surface brightnesses [(a) /io = 20 mag/arcsec 2 , (b) /j,q — 22mag/arcsec 2 , & (c) 
/i = 24 mag/arcsec 2 for exponential profile galaxies, and (d) (j,q = 20 mag/arcsec 2 , (e) fi = 22 mag/arcsec 2 , 
& (f) fio = 24 mag/arcsec 2 are for de Vaucouleurs profile galaxies], calculated using standard methods (heavy 
lines; eqn. |l]), and calculated using the proper corrections for lost light (light lines; eqn. [2]). The upper row is 
the luminosity function and V/V max distribution derived from the entire sample. The middle row is for the 
galaxies which fall below the mean redshift of the sample, and the lower row is for the galaxies falling above 
the mean redshift; the mean redshift changes with surface brightness. For the cases where the luminosity 
function appears to evolves with redshift, Schechter luminosity function fits are plotted superimposed to 
facilitate comparison; the dotted line in the lower left panel is the luminosity function fit to the low redshift 
data. All samples have 5,000 galaxies, and assume a limiting isophotal magnitude of mu m — 25 mag/arcsec 2 , 
measured within an outer isophote of \xu m = 22 mag/arcsec 2 , and 1" seeing, with a (3 — 5 Moffat profile 
distribution. They do not include the effects of random galaxy inclinations. Note that the luminosity 
function has a higher normalization for the lower surface brightness galaxies, due to the systematically 
lower accessible volume, and thus the smaller mean sample redshift (Figure 0; also Phillipps et al. 1990 and 
references therein). 

Fig. 4. — The first moment radius r\ measured above some limiting surface brightness fiu m (upper panel) 
and the resulting error in the absolute Kron magnitude measured within 2r\ (i.e. AM = —2.5 lg /(2ri), lower 
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panel), as a function of A^ = fj,u m — /Xo, for different angular scale lengths a = 0.25", 1", & 2" in the presence 
of a Moffat PSF with FWHM = 1", for [a] exponential disks and [b] de Vaucouleurs' profile ellipticals. Lower 
surface brightness galaxies (small A/x) have systematically underestimated Kron magnitudes, because the 
value of ri is underestimated when only a fraction of the galaxy's profile is used to calculate r±. Elliptical 
galaxies, which have a larger fraction of light at large radii, are subject to larger errors in the measured 
luminosity than exponential disks. The curves are purely theoretical, and do not take noise and other forms 
of measurement error into account. For reference, when measured within r — > oo, ri(oo)/a = 2 for a pure 
exponential profile, and ri(oo)/a = 3.8424 for a pure de Vaucouleurs' profile. The biases are identical for 
the same for galaxies with the same ratio of the FWHM to the scale length a. 

Fig. 5. — The error in the measured absolute magnitude as a function of redshift, for aperture and 
Petrosian magnitudes (top and bottom, respectively), for both exponential (left) and de Vaucouleurs 
(right) profile galaxies, in the presence of a 1" FWHM Moffat profile PSF with (3 = 5 (comparable to 
the middle row of Figure [[]). As in Figure 0, the different lines represent galaxies of different physical 
sizes (a = 0.25 kpc, 1 kpc, 3 kpc: dashed, dotted, and solid lines, respectively). Because both aperture and 
Petrosian magnitudes scale linearly with surface brightness, the given plots of AM are independent of galaxy 
surface brightness. Note the factor of 5 difference in scale between the aperture magnitude plots and the 
Petrosian magnitude plots. 
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Figure |[c] 
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